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I. INTRODUCTION 



In paper [1] we have shown that, by choosing the extensive variables (U, V, X 1 , . . . , X n ) as independent vari- 
ables, and by postulating that the integrablc Pfaffian form 5Q rev representing the infinitesimal heat exchanged 
reversibly is homogeneous of degree one and that the homogeneity symmetry is non-trivial [for the definition of 
non-trivial symmetry see the next section], it is possible to find immediately an integrating factor for SQ rev . In fact, 
if Y = U du + V dv + X)i X' 1 dx* is the Liouville vector field associated with the homogeneity symmetry, one has 
that an integrating factor is given by / = SQ rev (Y) = iySQrev ^ (see [1] and references therein). The entropy 
S(U, V, X 1 , . . . ,X n ), which represents the fundamental equation in the entropy representation, is then straightfor- 
wardly recovered. The framework presented in paper [1] can be generalized in the following sense. A generalized 
thermodynamics where SQ rev is a quasi-homogeneous integrable Pfaffian form is proposed. Quasi-homogeneity is 
realized to be a property which characterizes the behavior of thermodynamic functions like e.g. the entropy in the 
black hole case, in some statistical mechanical models involving non-relativistic matter interacting via the Newto- 
nian potential and also in the case of self-gravitating radiation. When gravity plays an important role the standard 
extensivity property of thermodynamics does not hold; nevertheless, one can still find a form of thermodynamics in 
which each thermodynamic variable follows a power scaling law where the power can be different from one or zero, 
i.e., the variables are allowed to be neither extensive nor intensive. Black hole thermodynamics is the most evident 
and special thermodynamics belonging to this framework. 

The plan of the paper is the following. In sect. II we give some definitions and then in sect. Ill, we propose a 
theoretical framework for quasi-homogeneous thermodynamics, by generalizing the standard thermodynamics case 
[1]. In sects. IV, V and VI examples displaying a quasi-homogeneous thermodynamics are discussed. In sect. IV, 
the role of quasi-homogeneity in black hole thermodynamics is shown; in sect. V, we discuss a model of fermionic 
non-relativistic matter with Newtonian interaction which displays a quasi-homogeneous behavior. In sect. VI, the 
general relativistic case of thermal geons and self-gravitating radiation are discussed. Sect. VII concerns the discussion 
of the physical consequences for thermodynamics when the requirement of homogeneity is relaxed. In sect. VIII a 
summary of constructive assumptions is made. In sect. IX the generalized Gibbs-Duhem identities are discussed. In 
sect. X the energy representation and the Legendre transform are analyzed. In sect. XI some further suggestions 
from black hole thermodynamics are discussed. Sect. XII has an heuristic nature, a link between the thermodynamic 
limit and the formalism of quasi-homogeneous thermodynamics is proposed. In the Appendix further mathematical 
and physical aspects are explored. 



II. QUASI-HOMOGENEOUS PFAFFIAN FORMS 

Given a set of real coordinates X - — (^X y . . . y X ) and a set of weights a = (a\, . . . , a n ) £ IR™, a real- valued function 
G(x 1 , . . . , x n ) is quasi-homogeneous of degree r and type a. if, under the action of the one-parameter group of quasi- 
homogeneous dilatations 

g T : (x\...,x n )^(e aiT x\...,e a " T x n ), (1) 

where r £ IR, one finds 

G(g T x) = e rT G(x). (2) 
A differentiable function G on a open connected domain 

V = {x £ V\g T x £ V V x £ V, t £ IR} C IR™ (3) 
is quasi-homogeneous of degree r if and only if [2,3] 

D G ~ r G, (4) 

where the Euler vector field 

D = a lX i A + . + 0lnX n A_ (5) 

is the infinitesimal generator of the transformation. Notice that a quasi-homogeneous transformation is also called 
"similarity transformation" and "stretching transformation" [4] . The identity (4) is a generalization of the Euler iden- 
tity for homogeneous functions. Quasi-homogeneity for functions and vector fields is defined in Ref. [2]. Notice that, 
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according to the definition given above, one cannot avoid specifying, together with the degree of quasi-homogeneity, 
the type of the quasi-homogeneous object one is considering. 

Homogeneous functions are a subset of quasi-homogeneous functions, in fact they have all the weights equal to one. 
We define strictly quasi-homogeneous functions the functions which satisfy the above definition with weights which 
cannot be all equal. A quasi-homogeneous function of n variables is characterized by n weights and its degree. If all 
the weights are undetermined and arbitrary, or if they are in part undetermined and in part equal to one, then all the 
weights can be put equal to one, and then the function is actually homogeneous. The viceversa is not true, because 
there are homogeneous functions which don't admit different weights (i.e., they are not strictly quasi-homogeneous). 
A trivial example is given by g(x,y) — x 2 + y 2 , which is homogeneous of degree two but it is not also strictly quasi- 
homogeneous. Instead, the function h(x, y) — xy is homogeneous of degree two and it is also quasi-homogeneous with 
undetermined weights. In the following, "quasi-homogeneous" is used as synonymous of strictly quasi-homogeneous, 
unless a more general discussion is required. 
A Pfaffian form 



E 



i(x) dx l (6) 



1 . ,x™ i-> A" 1 x 1 , . . . , A Q " x n (7) 



is defined to be quasi-homogeneous of degree r e IR if, under the scaling 

x 

one finds 

LO \ T LO. (8) 

This happens if and only if the degree of quasi-homogeneity deg(w i (x)) of u>i(x) is such that deg(wj(x)) = r—on Vi = 
1, . . . , n. Let us assume that a quasi- homogeneous u> is integrablc, i.e. cuAdui = 0. A quasi-homogeneous transformation 
is a symmetry for w (see e.g. [5,6]), in the sense that 

(L D u) A uj = (9) 

where is the Lie derivative associated with D defined in (5). An integrating factor / such that the form uj/f is 
closed can then be constructed by contracting the vector field D with lo: 

f = uj{D)=i D {w), (10) 

and w(D) ^ is to be assumed. This means that the symmetry generated by D is nontrivial, in the sense that it 
does not belong to the distribution of codimension one which is associated with the kernel of lo. Then, the symmetry 
is not tangent to each leaf of the foliation associated with the integrable 1-form u, but it carries leaves onto other 
leaves. In this sense, a nontrivial symmetry is transversal with respect to the foliation. 
One gets 

^)=0. (.1) 

The integrating factor / is a quasi-homogeneous function of degree r, because D «d(^) = in L^lu = r id^jj) [due to 
the Cartan formula Lx iy — iy Lx = i[x.y] an d to i = 0]. Then u/f is of degree zero and its integral can be found 
only by quadratures (cf. sect. Ill); one gets 

S-So^j r j- (12) 

It can be shown that there exists a quasi-homogeneous function F of degree one with respect to D such that 

S = log(F) (13) 
D F = F, (14) 

thus it holds DS = 1 and S is not quasi-homogeneous. On this topic see Appendix A, where the case of a generic 
transversal symmetry is treated. The above construction is completely analogous to the construction developed for 
homogeneous integrable Pfaffian forms in Rcf. [1]. 
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III. A GENERAL FRAMEWORK 



Quasi-homogeneous thermodynamics [recall that we mean strictly quasi-homogeneous] overcomes the standard 
distinction between extensive and intensive variables. We propose the following generalization for the intcgrablc 
Pfaffian form 5Q rev which represents the infinitesimal heat exchanged by the system: 

n 

SQ rev = dU*-J2 Ci dX l * (15) 

i=l 

where the asterisk indicates that both the independent variables U*, X 1 *, . . . , X n * and the dependent ones £*,...,£* 
in (15) are generalizations of the usual ones in the sense that they are not simply intensive and/or extensive but 
quasi-homogeneous, in such a way that the Pfaffian form SQ rev is quasi-homogeneous of degree r. Moreover, for the 
sake of definiteness and in order to fit some requirements for the validity of Frobenius' theorem, SQ rev is assumed to 
be of class at least C 1 in the thermodynamic domain (except, maybe, on the boundary as e.g. the surface T = 0). 
The thermodynamic domain D is assumed to be a simply connected set which satisfies (3). A further requirement 
for T> has to be introduced if the entropy S* is required to be superadditive: T> has to be closed with respect to the 
sum (see sect. VII for further details). Cf. also subsect. IIIB1. [In standard homogeneous thermodynamics, the 
thermodynamic domain can be assumed to be a convex cone [1], which is a cone with the property to be closed under 
addition ( [7], pp. 13-14). One can also assume that the domain V is still a convex cone with the further requirement 
that it has to be invariant under quasi- homogeneous transformations. This can be obtained by considering a set 
C C IR" +1 which is invariant under quasi- homogeneous transformations and then the set K, of all the positive linear 
combinations of elements of C. The set JC is the smallest convex cone containing C ( [7], p. 14). This cone JC is 
trivially still invariant under quasi-homogeneous transformations.] 

In the following, it is useful to refer to the variables U* , X 1 *,..., X n * as would-be extensive variables, and to the 
variables £jf , . . . , £* (and also T*) as would-be intensive variables. The variables U* , X 1 *, . . . , X n * could also be 
chosen to be such that they are additive, i.e., if one considers a system which is composed by two non- interacting 
subsystems, then X* = XI + X% . See also Ref. [8]. 

We assume that the X % * are of degree af, notice that the degree of U* is a = r, i.e., the degree r of SQ rev and the 
weight a of U* have to coincide in every case. The Euler vector field is 

D = aU* t^+Y oaX 1 * -4-. (16) 
dU* ^ dX % * v ' 

i 

An integrating factor for SQ re v is given by 

f* = aU*-Y,OL i QX i *, (17) 

i 

and it is assumed that /* ^ as in standard thermodynamics and, moreover, it is assumed that /* > 0, which is 
related to the positivity of the absolute temperature. Then, it is possible to show that, as in standard thermodynamics, 
the potential S*, which is not quasi- homogeneous, is associated with a positive definite potential S* which is a quasi- 
homogeneous function of degree one with respect to the Euler vector field D: 

*-* = / r # = "*(£)• (18) 

The proof is found in Appendix A. 

Before continuing our analysis, we recall that a detailed and important study on quasi-homogeneous functions and 
their application to scaling and universality is the subject of Ref. [9,10] [therein, quasi-homogeneous functions are 
called "generalized homogeneous functions" , which is a better denomination, but we adopt a common mathematical 
denomination]. See also [11]. Further mathematical properties are found in Ref. [12] [note that, therein, quasi- 
homogeneous functions are called "almost-homogeneous functions" ( [12], p. 231), whereas the definition "generalized 
homogeneous functions" is assigned to a further generalization of the equation defining homogeneous and almost- 
homogeneous functions ( [12], p. 304)]. 

A. metrical entropy 

By analogy with the construction for standard thermodynamics, a possibility is that S* is the metrical entropy for 
the system. This is what happens in standard homogeneous thermodynamics and also in black hole thermodynamics, 
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and we conjecture that this occurrence is not special. See however the discussion at the end of this section. 
We define T* > by means of 

dS* 1 

QU~, = > °> (19) 

which means that S* is assumed to be monotonically increasing in U*. T* is another integrating factor for (15), it is 
such that 

SQ rev = T* dS*. (20) 

T* is a quasi- homogeneous function of degree r — 1 [see also theorem 1 of [9], where it is shown that the partial 
derivative dg/dX of a quasi-homogeneous function g of degree a with respect to a variable X of weight b is a quasi- 
homogeneous function of degree a — b and the same type as g\. When r ^ 1, the function T is not intensive, but 
changes under quasi-homogeneous rescalings of the system. In the case of a Kerr-Newman black hole, one finds e.g. 
that, by doubling the mass and the charge, and by quadruplicating the angular momentum, the temperature becomes 
one half the temperature of the original black hole state. This behavior is well-justified in the light of the Hawking 
effect. 

As a consequence of (19) and of (18), one finds 



dS* 1 1 dS* 



dU* f* S* dU* T*S* 

thus 



(21) 



f*=T*S*, (22) 
as in standard thermodynamics. From (20) it is easy to show that the quasi-homogeneous entropy can be written as 

i 

In fact, T* is an integrating factor of degree r — 1 for 5Q rev , and dS* = SQ rev /T* is an exact quasi-homogeneous 
Pfaffian form of degree one. Then, as a consequence of lemma 1 in Appendix A (cf. also result 1 therein), one finds 
S* = io (SQrev/T*), i.e. (23) holds. When all the weights in (23) are equal to one, the well-known expression for the 
homogeneous thermodynamic entropy is recovered. 



B. thermodynamic foliation 

As in the case of standard thermodynamics, one can require that the thermodynamic foliation is defined by the 
leaves S* — const, everywhere in the thermodynamic manifold, except maybe on the boundary. Singularities for S* 
can occur where /* = 0, i.e., in the set Z(f*) = Z(T*) U Z(S*). 

The surface T* = is expected to represent an adiabatic boundary of the thermodynamic manifold. Notice that, as 
in the case of standard thermodynamics, the set Z(S*) of zeroes of the function S* , if non-empty, is assumed to be 
contained in the set Z(T*). In standard thermodynamics the occurrence of a zero for S at a temperature T > can 
be considered pathological, a system in such a state should necessarily absorb heat in a neighborhood of this state, 
whichever reversible process could be considered. Cf. [1]. The same would happen in the quasi-homogeneous case. 

Because of (18), the singular values for the thermodynamic foliation are represented by the points where S* = (if 
any). 



1. a stronger assumption on the domain 

One could also assume that 

/* > U* > b(X x *,..., X n *), (24) 

where b(X 1 *,..., X n *) is a quasi-homogeneous function of degree a and weights (ct\, . . . , a n ). This function 
b(X 1 *,..., X n *) plays the role of lowest energy for the system, because by definition U* > b. The domain V 
has to include the following set: 
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epi(6) ee {([/*, X 1 *,..., X n *)\{X 1 *,..., X n *) e /C ( „), U* > ft^ 1 *,..., X" *)}, 



(25) 



where /C(„) C IR" is an open connected set. epi(&) is by definition the epigraph of the function b. If /C(„) is closed 
under quasi-homogeneous dilatations (X 1 *,..., X n *) i— > (A Ql X 1 *, . . . , A Q ™ X™ *), then epi(fe) is closed under quasi- 
homogeneous dilatations involving also U*, because 6(A Ql X 1 *,..., A a ™ X"*) = A" b(X x *,..., X n *) < A" [/*. 
Moreover, if /C(„) is also closed under the sum and if b is subadditive, then epi(6) is closed under the sum too, thus the 
closure of /C(„) under quasi-homogeneous dilatations and under the sum allows us to choose V = epi(6), in analogy 
with the standard homogeneous case. [Notice that if e.g. JCr n \ = IR™ , then it is a convex cone which is invariant under 
quasi-homogeneous dilatations] . An example where the domain V coincides with the epigraph of a quasi- homogeneous 
function is furnished by black hole thermodynamics. 

Under this assumption on the structure of the domain, one finds that it is not possible to find a state having S* = 
at T* > 0, because by introducing 

B* = U* - biX 1 *,.. .,X n *) > (26) 

one finds 

S*(B*,X 1 *,...,X n *) = S*(0,X 1 *,...,X n *) + dY T , {Y>xl } ^ >Xn ,y ( 27 ) 

which cannot be zero for any B* > (in the last formula mathematical conditions ensuring a continuous entropy at 
T* = have been implicitly assumed. Cf. [1]). 

As far as the reference state is concerned, we assume that states like the one where all the variables U* , X 1 * , . . . , X n * 
are zero, and any states which imply the absence of the system, are unphysical. The thermodynamic description starts 
being meaningful if a statistically relevant number of degrees of freedom is available. Cf. sect. II of paper [1]. 



C. possible ambiguities 

One may wonder what should happen if one finds that w is quasi-homogeneous but the exact weights are not a 
priori known. In particular, one can assume that their ratio is known, i.e., a/ at and aj/ai are known for all i ^ j; 
then, there is an overall undetermined multiplicative constant q ^ 0, in the sense that these ratios don't change if 
one multiplies all the weights by the same constant. One could e.g. know = aja\ and r^ 1 = aj/ai for a\ ^ 

and for all j ^ 1, but the absolute weight a\ is not known, thus the weights are (ai r^, ai, ai ...,ot\ r^). 
An example where this ambiguity appears is represented by black hole thermodynamics (see sect. IV). One can also 
assign the weights with respect to a scale A whose absolute weight could be unknown. The absolute weights are then 
of the form (q ot,q ai, ... ,q a n ) with q undetermined. We treat the latter case without loss of generality, being the 
former equivalent to it under the map ai i-»gai, q u\ q a, q ol\ r'p q ai for all i ^ 1. As a consequence of 

the undetermined overall factor q, one has a one-parameter family of Euler vector fields {D^}, with 

D M=^U*-^ F +J2^iX i *^-. (28) 

i 

Thus, there exist a one-parameter family of integrating factors, with f*^ = SQ rev (D( q ^), and a one-parameter 

family {S( q ^} of potentials such that 

SQ r 



Crev 



and a one-parameter family of entropies {-S^)} satisfying 



dS {q) (29) 



dS (a) 

D («)%=%- ( 31 ) 



G 



The true metrical entropy 5** belongs to the family {S^} but is undetermined because the overall factor q is unde- 
termined. In order to see which relation occurs between the various entropies in this family let us consider a fixed 
value q of the parameter q. Then we obtain the entropy S^ q y If one considers a generic q ^ q, one has 

%) = \ ( 32 ) 

and 

A« = |/(V ( 33 ) 



Then, 

SQrev q SQ 

rci 

f*q) 9 f(q) 



(34) 



diog(sy = Uio g (sy. (35) 



As a consequence, by assuming in the integral of SQ rev / f the same reference state (Uq , Xq *,..., Xfi *), one finds 

(«) ^ ' '•••' A o J 

(skills, xs*)W 



S^U^X 1 X n *)= ^f ( ^° y f" '"'"'^L (SfaiU^X 1 A" *)f/\ (36) 



i.e. 



where Cg,g = Sf-JU^ ,Xq *,..., Xft *) / (ST-JUq , Xq *,..., Xfi *))«/« is a constant, once one has fixed the reference 
state. The metrical entropy S* is related with the entropy S*-^ by a simple power law. Cf. also Appendix A. 
The one-parameter family of temperatures {T* q ^}, with 

1 dST n) 

- {q) (38) 



is such that 



^- = Cw| W 75 " 1 ^- ( 39 ) 

Analogous relations exist for the other partial derivatives of and Sf-y Notice that it holds f* q ^ = S* q y 
because 

^k-^-^^k-^ (40) 

W f* q) S* {q) dU* T* q) S* q) - ^ 1 



1. phenomenological resolution of the ambiguity 

The metrical entropy S* could be phcnomcnologically identified by means of an absolute temperature thermometer, 
made of standard matter, allowing to find, as the parameters U* , X 1 *,..., X n * are varied, the function T* describing 
the absolute temperature of the system. As a consequence, also S* can be recovered. We remark that the resolution 
of the above ambiguity by means of a phenomenological input is completely in the spirit of thermodynamics, which 
has a phenomenological nature as far as statistical mechanics is not taken into account. 
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2. how to fix the absolute weights by means of statistical mechanics 



A priori, even in statistical mechanics, there is a scaling ambiguity, unless the absolute weights are somehow 
fixed. On this topic, see also sect. XII. In fact, let us assume that the statistical mechanical entropy [calculated 
by means of some sort of thermodynamic limit or by means of some mean field approximation; cf. sect. XII] under 
rescaling by means of a scale A is a quasi-homogeneous function of degree as and weights a\, . . . , a n . This means 
that S sm (X a U, X ai _X 1 , . . . , A"" X n ) = X as S sm (U, X 1 ,..., X n ). By changing the scale A = A 9 , the scaling of the 
variables becomes X q ai . Thus, there is an overall factor q of ambiguity. Such a redefinition of the scale, unless 
the absolute weights are somehow fixed, allows to re-map any quasi-homogeneous function of degree a^O into a 
quasi-homogeneous function of degree one, as noticed in [10]. The statistical mechanical entropy is then assumed to 
satisfy 

5 sm (A« Q C/,A« Q1 X\...,X" a " X n ) = X qas S sm (U,X\...,X n ). (41) 

where the overall factor ambiguity is enhanced. With respect to the Euler vector field which corresponds to the 
above weights it holds D( 9 ) S sm — q as S sm . 

We know that from a thermodynamic point of view there is a one-parameter family of possible metrical entropies, each 
of which is of degree one. The statistical mechanical entropy S sm has to coincide with the metrical thermodynamic 
entropy S* and it has to hold T dS sm = T* dS* — SQ rev . As a consequence, the statistical mechanical entropy has to 
be a degree one quasi- homogeneous function. This requirement fixes q as = 1, and fixes unambiguously the weights 
[then, it eliminates also the ambiguity in the one-parameter family of thermodynamic entropies]. 

We could also relax our conjecture about the identification of S* in formula (18) as the metrical entropy The 
Clausius-like formula (18) gives in general an empirical entropy which is to be uniquely related with the metrical 
entropy by identifying the absolute temperature scale as in the standard procedure. One could still assume that the 
metrical entropy is quasi-homogeneous, but it would be possible in line of principle to require that its degree is q, with 
q non-necessarily equal to one. Then, even in the case where all the weights are known, one would find 5^ ctric cx (S*) g , 
where S* is obtained from formula (18). Cf. also Appendix A. The above phenomenological procedure could be still 
allowed. Cf. also sect. XII and the conclusions herein. 



A further consideration is in order, and it concerns the possible presence of variables of weight zero. All the math- 
ematical approach developed in this section remains unaltered. It is to be noted that, if a variable I* having weight 
zero is present, the term I* d/I* does not appear in the Euler operator and in the integrating factor / the term 
proportional to I* is missing. 

Notice that it is possible that the thermodynamic description allows a reducing/enlarging of the thermodynamic 
space as the one discussed in Ref. [1]. This means that part of the variables X % * appearing in 5Q rev and in S* could 
be put equal to zero consistently, which means that a meaningful thermodynamic description is still allowed when 
X % * = [in standard thermodynamics such variables play e.g. the role of external fields [1]]; Consistency requires 
that, if X 1 * = is allowed, then £* = as X % * — [1]. Black hole thermodynamics represents a good example, and 
this topic is discussed in sect. XL 

We underline that this kind of generalization of thermodynamics can be considered as a sub-case of a generalization 
proposed by P.T.Landsberg in Ref. [13]. Therein, a discussion of generalized ideal gases appears, where the properties 
of the following generalized entropy are analyzed: 

UV 9 

S = bN\og(a 1 ^ r ). (42) 

Given the Euler field 

D^-cvJL+PV^+inJL, (43) 

it is easy to show that S in (42) is quasi-homogeneous (of degree 7) if and only if a + g (3 — /i 7 = If the latter 
constraint is not implemented, then S is not quasi- homogeneous. Homogeneity requires h = g + 1. 
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IV. THE BLACK HOLE CASE 



We summarize herein the results obtained in Ref. [14] in the black hole case (we put h = c = G = kb = l, where kb 
is the Boltzmann constant; moreover, we work with unrationalizcd electrical units). For a general discussion about 
black hole thermodynamics see e.g. [15,16]. Herein, as well as in sect. V and in sect. VI, we do not use the asterisk, 
which is in general used only in theoretical sections. The Pfaffian form of black hole thermodynamics for black holes 
of the Kerr-Newman family is 

SQ rev =dM -QdQ-tt dJ, (44) 

where the parameters M, Q, J are the mass M, the angular momentum J and the charge Q of the black hole; the 
angular velocity 

n = T7 ~, ( 45 ) 

M 2 M 2 — Q 2 + 2 M V / M 2 -Q 2 -J 2 /M 2 

and the electric potential 



* = Q(M+VM 2 -Q 2 -J 2 /M 2 ) 

2 M 2 - Q 2 + 2 M y/M 2 -Q 2 - J 2 /M 2 

are associated with the black hole horizon. M in (44) plays the role of internal energy, the remaining terms in (44) 
being "standard" work terms. Each black hole of the Kerr-Newman family is a stationary solution of the Einstein 
equations, and it can be considered, from the thermodynamic point of view we develop here, rather loosely speaking, 
as an "equilibrium state of the geometry". The variation occurring in (44) is taken along such stationary solutions. 
The thermodynamic domain for the non-extremal black holes is the (open) non-extremal manifold implicitly defined 
by the inequality M 4 - M 2 Q 2 - J 2 > 0; the extremal sub-manifold M 4 - M 2 Q 2 - J 2 = is a boundary of the 
former, and its status is analyzed in Ref. [17] (therein a discussion of the third law of thermodynamics also in the 
light of this approach is found. We limit ourselves to recall that in black hole thermodynamics the entropic form of 
the third law is violated [18] but still the unattainability of T = holds under suitable hypotheses). 
The Pfaffian form (44) is integrable, it satisfies the condition SQ rev A d(SQ rev ) = 0, i.e. 

-dj $ + d Q Q + $ d M - n d M $ = 0. (47) 

Moreover, the Pfaffian form (44) is quasi-homogeneous. In fact, under M — > A"M; Q — > X a Q; J — > A 2 " J, one obtains 
SQrev — * X a 5Q r ev, where A > 0. The weights a, a, 2a are not known, only their ratio is determined. Then, the degree 
of SQ rev is not fixed [14]. Nevertheless, a thermodynamic construction is allowed. One finds that the integrating 
factor for 6Q re v is given by / = a (M — $ Q — 2 ft J) , that is 

/ = a VM 2 -Q 2 - J 2 /M 2 ; (48) 

[it is interesting to notice that this integrating factor is proportional to the horizon coordinate c introduced by B.Carter 
in [19]. Thus, it is constant on the horizon]. Then, one gets [14] 



/ 



^"^5-5o = i-logfA), 



/ u 2a °\A 



io 



where A = 4 tt (M 2 (1 + y/l - Q 2 /M 2 - J 2 /M 4 ) 2 + J 2 /M 2 ) is the black hole area and is positive definite. The 
main difficulty is in finding the metrical entropy, because of the undetermined overall multiplicative constant a. 
Both a = 1/2 and a undetermined proportionality constant Co appearing in the Bekenstein-Hawking law S — cq A 
are recovered e.g. by comparing the temperature of the black hole with the Hawking temperature, and one finds 
c = 1/4; also the Euclidean path integral method can be used with this aim [20] (see also the papers collected in 
[21]), as well as a phenomenological plot of the temperature as a function of M, Q, J. [Moreover, we notice that the 
same result S cx A could be obtained by requiring that the metrical entropy is quasi-homogeneous of degree one with 
respect to the Euler field D which is obtained by normalizing to one the biggest weight 2a: 

1 d 1 d d 
D =2 M OM + 2 Q dQ +J OJ- ^ 
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But this normalization rule should be justified on a general thermodynamic footing, and the case J — (Reissner- 
Nordstrom solution), being homogeneous in M, Q, shows that the physically correct value a = 1/2 cannot be found 
on the ground of such a rule] . 

S is quasi-homogeneous of degree one and of type (1/2, 1/2, 1) and it is well defined everywhere in the thermodynamic 
manifold. It is such that 5 = only for M = Q = J = 0, which should not be considered as a state belonging to 
the thermodynamic manifold [see the discussion in section III; moreover, near M = quantum gravity effects are 
non- negligible and General Relativity is expected to be non-viable]. Thus, Z(S) = and one finds Z(f) = Z(T), and 
S is regular everywhere; as a consequence, S and S could be continued even on the extremal submanifold (but see 
[17] for arguments against this continuation). All non-extremal states have T > 0, whereas the extremal boundary 
corresponds to T = 0. As it is well-known, one obtains an entropy which is not concave, but simply strictly superad- 
ditive when the merging of two black holes is considered. This ensures the second law of thermodynamics [13,22,23]. 
See also sect. VII. 

By introducing 



which is a quasi-homogeneous function of degree 1/2 and weights (1/2, 1), as it can be easily verified (see also [17]), 
the domain can be also written as follows: 



In fact, M — b(Q, J) is the physical [i.e. having positive mass] solution of M 4 — M 2 Q 2 — J 2 = and the black hole 
manifold is M A — M 2 Q 2 — J 2 > 0, which means M > b. The function b(Q, J) indicates, for given values of Q, J, 
which is the lowest mass such that a black hole solution can exist. The lowest mass coincides with the mass of the 
extremal black hole having Q, J as charge and angular momentum respectively. 

We point out that the quasi-homogeneous behavior of black hole thermodynamics is not a special feature of the 
Kerr-Newman family. In fact, it can be realized easily e.g. in the case of the Kerr-Newman-Anti-De-Sitter (KN-AdS) 
case [24]. A peculiar feature arises, one has to define a thermodynamic angular velocity as the difference between 
the angular velocity of the horizon and the angular velocity at infinity, and the same is true in the case of the electric 
potential [24]. Nevertheless, these definitions are necessary, e.g. the "electric" work term has to involve the difference 
between the potential at the horizon and the potential at infinity when the latter is not zero, and, on the same footing, 
the redefinition of f2 is due to the fact that the angular velocity is not zero at infinity. The change in the energy 
associated with these work terms has to involve such a difference, which does not appear in the KN case because of 
the vanishing of $ and of f2 at infinity. 

General Relativity is surprising from a thermodynamic point of view. Black hole solutions of the Einstein equations 
are involved with an integrability condition (47), and they allow a fine explicit thermodynamic construction a la 
Caratheodory. One finds a link with thermodynamics which is a priori unexpected and it is not simply the formal 
analogy between laws of thermodynamics and laws of black hole mechanics, because the above construction belongs 
to the thermodynamic framework, apart from the geometric inputs f2, which come from General Relativity, and a 
comparison with the Hawking effect, which fixes a = 1/2. See also the discussion in Ref. [14]. 

V. NEWTONIAN GRAVITY. A MODEL OF HERTEL, NARNHOFER AND THIRRING 

There exists an interesting statistical mechanical calculation for gravitating fermions which corroborates the idea 
that, in presence of gravity, a quasi- homogeneous thermodynamics is allowed [25-27]. The model (HNT model 
in the following) involves non-relativistic fermions interacting by means of a Coulomb potential and a Newtonian 
potential, with vanishing [or small] total charge. Thus, the calculation does not involve General Relativity but only 
Newtonian gravity. Nevertheless, the result is still interesting and, to some extent, puzzling. The quantum mechanical 
Hamiltonian which is considered is 




(51) 



V = {(M, Q, J)\(Q, J) e IR 2 , M >6(Q,J)} = epi(6). 



(52) 




(53) 
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The authors state that a thermodynamic limit exists if the following asymptotic behavior is allowed for U, V: V <~ 1 /N, 
U ~ N 7 / 3 , and the micro-canonical entropy is of degree one in the number TV of fermions. The scaling properties of 
U, V with N are peculiar, they are homogeneous of degree one if the potential term in (53) is suppressed by means 
of a multiplicative factor 1/N 2 / 3 , which means that, as the system gets larger, the interaction becomes weaker [27]. 
If this suppression is not in agreement with the physics, and if the fermions do not become relativistic, then the 
micro-canonical entropy S(U, V, N) is a quasi-homogeneous function of degree one, and the weights of U, V, N arc 
7/3, —1, 1 respectively. Cf. also sect. XII, where the notion of quasi-homogeneous thermodynamic limit is introduced. 
The Euler field is 

DuvN = lujL- V JL +N jL. (54) 

For details of the highly non-trivial calculations, see Refs. [25,26,28]. We are interested here in pointing out that, 
as stressed in Ref. [25], there is a one-parameter family of equivalent limits N — > oo, whose parameter 7 remains 
fixed by the requirement that the ground state energy goes like N 7 / 3 , according to a result of Ref. [29]. A scaling 
ambiguity, which would reflect itself in the undetermined weights (—27 + 5/3, 37, 1) of (U, V, N) respectively, is thus 
solved by comparing with the scaling behavior of the ground state. See also [28]. Another important point is that 
the temperature-dependent Thomas-Fermi equation, which allows these scaling properties, becomes exact in the 
thermodynamic limit for the system under discussion [25]. 

We limit ourselves to underline that Newtonian gravity furnishes a statistical mechanical example where unusual 
scaling laws for the "extensive" variables are allowed. The purely attractive nature of gravity plays a major role, 
because it does not saturate, i.e., it does not allow to obtain a ground state proportional to the number of particles 
[29]. Notice also that a lack of concavity is allowed [27]. 

It is interesting to notice that another model involving different scaling relations has been developed in Ref. [30]. 
Therein, a classical gas of non-relativistic particles which interact by means of Newtonian gravity is considered in a 
diluted regime where particles of mass m are enclosed in a box; the following behavior is recovered in a non-conventional 
thermodynamic limit where V/N 3 -^const. as N — > 00. A dilution parameter 

- 1 UV "' (55) 



G m 2 N 2 

is kept constant in the thermodynamic limit. One finds that S(U, V, N) is quasi-homogeneous of degree one and 
weights (1,3,1), as can be inferred from Ref. [30]. Cf. also sect. XII. This model differs from the previous one 
because it involves classical matter, the gas is kept diluted in the thermodynamic limit, no collapse is considered. 
Note that even the HNT model can be related with the parameter £ introduced in (55), in fact even in HNT model 
this parameter is constant in the appropriate thermodynamic limit. £ can be qualitatively interpreted as the ratio 
between the internal energy and the Newtonian gravitational energy, in fact the gravitational energy of a homogeneous 
sphere of mass M = N 2 m is proportional to GM 2 /V 1 / 3 . Thus, in both cases, the requirement of a non-negligible 
contribution of the gravitational energy to the internal energy is realized in the thermodynamic limit. In the HNT 
model, a self-bounded system is considered, like a (cold) star, where V contracts as the mass is increased; in the 
model of Ref. [30], the physical conditions arc different, a gaseous system is considered, the energy per particle 
remains constant and small, and the systems remains gaseous in the thermodynamic limit. 

In the regime where gravitational forces cannot be neglected, one can expect that a different kind of quasi-homogeneous 
thermodynamics can arise, with different weights, depending on the actual physical conditions, e.g., the amount and 
the type of matter. Semi- relativistic matter is expected to behave in a different manner (cf. Ref. [29]). The point 
is that a different scaling behavior of the ground-state of the matter can be expected if relativistic effects in the 
kinetic energy are non-negligible. A self-consistent calculation taking into account General Relativity is expected 
to give another kind of contribution. An approach, combining in a self-consistent way general-relativistic equations 
and finite-temperature Thomas-Fermi approximation for the thermodynamic functional would be required. In this 
direction, see Ref. [31]. It is also interesting to point out that such an approach probably would not be useful for 
understanding black hole thermodynamics, in fact no stable ground state for the matter could exist by hypothesis 
(a collapse of the matter beyond the possible stable states represented by white dwarf and neutron stars should take 
place, with no possibility to prevent the formation of an event horizon). 



VI. A LOOK AT THE THERMAL GEON AND AT SELF-GRAVITATING RADIATION 

We discuss naively the scaling properties of the so-called thermal geons [32] . A gcon is a "gravitational electromag- 
netic entity" which consists in a self-gravitating ensemble of electromagnetic modes [33] . The most stable configuration 
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is such that the photons are distributed within a toroidal region of the space-time, but also a spherical distribution 
can be allowed. Geons don't represent strictly stationary solutions of the Einstein equations, because a leakage of 
photons to infinity is allowed. Nevertheless, they can be considered as metastable quasi-stationary solutions. More- 
over, sizes are considered where no contribution of electron physics, due to pair creation by the electromagnetic field, 
is considered, and no zero-point energy is taken into account. In Ref. [32] the idea of a geon has been generalized 
to the case of thermally distributed electromagnetic modes. The main assumptions are the following ones: a) the 
metric is a spherosymmetric static diagonal metric in the Schwarzschild gauge; b) there exist two separate classes of 
electromagnetic modes; the first class is constituted by bounded null geodesies which represent modes whose energy 
cannot escape to infinity [the actual rate of flux to infinity goes to zero exponentially with the ratio between the 
dimensions of the geon and the wavelength A]. The second class is constituted by free electromagnetic modes, i.e. 
by null geodesies which reach infinity, c) The energy of the free modes is zero, the energy of the bounded modes of 
frequency f2 is [in natural units] 

which is the usual distribution for black body radiation and T is the temperature. See also [32] on this topic. The 
active region of the geon is defined as the region where the metric coefficient 500 satisfies |goo| < 1- By defining R as 
the radius of the active region and M as its mass, the following scaling laws can be deduced [32]: 

_ _ /15fiV\ 1/2 1 1 

*~^ s (wJ (57) 

As an order of magnitude estimation on the same footing of the previous evaluations, one can see that the entropy 
behaves as S oc VT 3 , where V is the volume of the active region, thus S cx T~ 3 . From the above naive considerations, 
one finds that S is a quasi-homogeneous function 

S(\- 2 M, \~ 6 V) = A~ 3 S(M, V); (59) 

In the case of a massless gas, there is an ambiguity in the identification of the absolute weights of the variables. There 
is an overall multiplicative factor to be determined in the Euler vector field. Let us assume that S is of degree one, 
as suggested by the formal picture of thermodynamics (cf. sect. Ill and sect. XII); then T has degree —1/3 and the 
scaling law for the fundamental equation is 

S(X 2/3 M, X 2 V) = \ S{M, V), (60) 

i.e., M has weight 2/3 and V has weight 2. We note that the ratio M/V 1 / 3 is scale-invariant and it is the same ratio 
which is kept constant in the thermodynamic limit discussed in [30]. 

The above result about thermal radiation can be confirmed by comparing it with the study of Ref. [34] . Black body 
self-gravitating radiation enclosed in a spherical box of radius R is considered in [34] . It is shown that the maximization 
of the entropy in a spherosymmetric geometry leads to the following results: i) the metric is static; ii) the perfect fluid 
of photons satisfies the Tolman-Oppenheimer-Volkov (TOV) equation. In particular, the TOV equation one obtains 
is scale- invariant. By introducing the density p(r), the mass m(r) = J Q r dyAiry 2 p{y) (r is the radial coordinate) one 
can define, by following [34], the dimensionless variables 

p = m(r)/r; (61) 
q = dm/dr = Anr 2 p(r); (62) 

and the TOV becomes equivalent to the following couple of equations: 

!«»• (63) 

No equilibrium is attained if p > 0.25 [34]. We are interested here in the scaling properties of the Pfaffian form 
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SQ rev = dM+p dV, (65) 

where M = m(R) is the mass, V is the volume (as seen from infinity) and p = p(R)/3 is the pressure as a function 
of the box radius. These variations are to be intended, as in the black hole case, as "on shell" , i.e., along static 
sphcrosymmetric equilibrium solutions of the TOV. Each solution represents a thermodynamic equilibrium state. 
Under the scaling M — > A M; V — > A 3 V; p — ► A~ 2 p, which corresponds to (30) of [34] and is allowed by the scale 
invariance of the TOV, one finds SQ rev — > A SQ rev . Particularly, one has 

D = a < M mi +3V w)> (66) 

where a has to be determined. From dV = AirR 2 dR and from (64) one finds 

SQrev = ^qdR (67) 



and the integrating factor 



SQ rev (D) =a(M + 3pV) = a ^ (3m + q) R. (68) 



Thus, 

u> 4 q dR 



f a 3/i + q R 



(69) 



By comparing this result with the ratio dS/S which can be obtained by the exact result (34) of [34], one finds 
that a = 2/3. As the consequence, the weights of D are the same as in the case of thermal geons (notice that 
the requirement for 5 to be quasi-homogeneous of degree one in the scale factor A would lead to the scaling r — > 
A 2/ 3 r . m ^ A 2/ 3 772* p — > A" 4 / 3 p in (30) of [34]). 



VII. CONSEQUENCES OF THE NON-EXTENSIVITY OF S 

We have to underline that the lack of extensivity has important consequences on the properties of the entropy In 
fact, if one considers for a continuous entropy 5 superadditivity (S), concavity (C) and homogeneity (H), one finds: 

a) (H) and (S) imply (C) [35]; 

b) (H) and (C) impiy (S) [36]; 

c) (S), with the condition S(0) = 0, and (C) imply (H) [8,22]. 

We recall that a function h(x) of n variables, collectively indicated with x in this section, is superadditive if 
h(x + y) > h(x) + h(y) for all x,y in its domain. Superadditivity in the case of entropy, means that the princi- 
ple of entropy increase in a irreversible adiabatic process is respected. Moreover, (C) requires a convex domain, (S) 
a domain which has to be closed under addition, and (H) a domain which has to be closed under multiplication of x 
by a real scalar A > (a cone). 

Notice that the condition 5(0) = appearing in c) is a mathematical condition which could be required even if 
the state x = does not belong to the thermodynamic space, in the light of our discussion in sect. Ill; see however 
also [22], where x = is included in the thermodynamic space. One could consider an extended domain including 
x = 0, on a purely mathematical ground, even if the thermodynamic formalism at such a point is meaningless. 
Moreover, 5 is required to be a continuous function. Then, if homogeneity does not hold and 5(0) = 0, either super- 
additivity or concavity has to be violated. It can be easily shown that the condition 5(0) = in c) can be replaced 
by the (natural) requirement that S(x) > for any state x in the thermodynamic domain. In fact, (C) implies 

S{\y+\z)>\s{y) + \s{z) (70) 

for any y, z; moreover, (S) implies 

5(2 x) > 2 S(x). (71) 
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If z = and y = 2 x, from (70) one obtains 

S(x)>ls(2x) + ±S(0). (72) 

Then, 1/2 5(2 x) + 1/2 5(0) < 5(x) < 1/2 5(2 x), which implies 5(0) < 0. As a consequence, the requirement of 
a nonnegative entropy 5, together with (S) and (C), implies 5(0) = 0, thus also (H) has to hold. This shows that 
the condition 5(0) = 0, in the framework of thermodynamics, where 5 > 0, cannot be actually considered as a true 
restriction leaving room for a thermodynamics in which a concave and superadditive but non-homogeneous entropy 
is allowed if 5 is defined in x = [of course, this holds as far as negative values of 5 are forbidden]. Notice that 
5*(0) = for a continuous quasi-homogeneous entropy defined in 0. 

Supcradditivity of the entropy means that for the thermodynamic system the entropy does not increase by frag- 
menting the system [8]. From the point of view of the energy representation, to be discussed in sect. X, superadditivity 
of the entropy is equivalent to the subadditivity of the energy under the condition (19) [for a proof, see Refs. [35,36]. 
This equivalence holds even if conditions (C), (H) are violated]. Subadditivity of the energy means that exploding 
is not energetically advantageous for the system [37,38]. The relevance of (S) for the second law of thermodynamics 
is discussed and underlined in Ref. [13]. Accordingly, one should privilege the superadditivity property against the 
concavity, and superadditivity should be required as a fundamental property for the quasi-homogeneous picture of 
thermodynamics. Notice that a lack of concavity for thermodynamics in presence of gravity is verified in Refs. [39,40], 
where negative heat capacities in presence of gravity are calculated. Therein, a discussion on how to deal with the 
lack of the standard stability properties of thermodynamics is sketched. See also [41,42]. It is to be noted that 
this lack of concavity in presence of gravity forces to abandon the usual Gibbsian scheme for thermodynamics, and 
the homogeneity property has to be withdrawn because of a), if one adopts the superadditivity as a fundamental 
property for the metrical entropy of a system. Due to the above equivalence between superadditivity of the entropy 
and subadditivity of the internal energy, one can naively justify this choice by underlining that, because of the purely 
attractive nature of gravity, internal energy of a self-gravitating system should be strictly subadditive, in fact energy 
decreases as the accretion of matter increases. Thus, there is an implosive instability of gravity, which eventually leads 
to the formation of black holes, to be identified as very long-lived metastable states. By passing, we note that there 
is another instability which is opposite in character with respect to the one implied by gravity. It involves system 
fragmentation/explosion because of a superadditive energy [37,38]. It should be characterized only by unstable states 
[systems with a suitable charge excess would be explosive [28]. See also the works of E.H.Lieb collected in Ref. [43], 
in particular the review [44]; see also [45]]. 

VIII. SUMMARY OF CONSTRUCTIVE ASSUMPTIONS 

We summarize the set of constructive assumptions upon which our approach is based. They correspond to the 
assumptions leading to standard homogeneous thermodynamics [1], the difference consists in the substitution of the 
homogeneity symmetry with the quasi-homogeneity symmetry and also in the explicit request that the entropy is 
superadditive. We comment only some assumptions which require further remarks with respect to the discussion 
developed in the text. 

al) The quasi-homogeneity symmetry which characterizes the thermodynamic system is translated into the quasi- 
homogeneity of the Pfaffian form SQ rev = dll* - Y%=i C dX l * . 

a2) The quasi-homogeneity symmetry of SQ rev is nontrivial. 

a3) The thermodynamic foliation is defined by the leaves 5* _1 (c), where c € IR is a constant, everywhere in V. 

This assumption does not leave room for foliations which are based on a quasi-homogeneous 5* which is allowed 
to be negative. Cf. the discussion in [1] for the homogeneous case. 

a4) The integrating factor /* is non-negative. 

a5) The metrical entropy 5* is quasi-homogeneous of degree one. 

This assumption could also be weakened, as discussed in the text. 
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a6) The metrical entropy S* is superadd itive. 

This assumption substitutes the requirement for a concave entropy occurring in standard homogeneous thermo- 
dynamics. Notice that in standard thermodynamics it would lead again to a concave metrical entropy, because of the 
implication a) of sect. VII. Cf. [35]. 

a7) We require that Z(S*) C Z(T*). 

This assumption is e.g. implemented as in sect. IIIB 1. 

One can also add the following assumptions which appear to be physically appealing: 
a8) We require that to each level set S* — const, corresponds a unique leaf. 

This assumption means that two states of the same quasi-homogeneous system lying on the same isoentropic surface 
are path-connected, i.e. it is possible to find a reversible adiabatic path connecting each other. 

a9) We require that J|tt is positive if V* is the volume of the system. 

This simply means that the pressure p* is positive whenever it is definite. 

alO) We require that /* = corresponds to an integral manifold of SQ rev . 

This assumption is natural and is related to the problem of the third law. Cf. [1] for standard thermodynamics; 
see also [17] for the black hole case and also Appendix B. 
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IX. GENERALIZED GIBBS-DUHEM EQUATION 



Quasi-homogeneity of S* allows to find a generalization of the standard Gibbs-Duhcm equation [46] . 



A. case where none independent variable of weight zero appears 

We first discuss the case where none independent variable of weight zero appears. By differentiating (23) one finds 
dS* = a -L dU* + a U* d (±) ]T a, £ dX* * - £ a t X" d (73) 
and by comparing with 

dS* = ± z dU*-Y J 7^dX i * (74) 

i 

one finds that the following generalized Gibbs-Duhem equation has to hold: 



This equation, as in the usual thermodynamic case, allows to express the differential of a would-be intensive variable 
as a function of all the other ones, and then to find such a would-be intensive variable apart from a integration 
constant (cf. [46] for the case of standard thermodynamics). In fact, one could e.g. find for 1/T*: 

rfl °S ( ^) = ~~~j7~~ dU* + ± ]T a, X* d(x^ . (76) 



Notice that 



d!og(^)=-dIog(r) + ^. (77) 

Eqn. (76) is implemented if (18) and (74) hold. Let us consider the inverse problem where one assigns n would- 
be intensive functions £*,...,£* which are quasi-homogeneous in U* , X 1 * , . . . ,X n * and such that SQ rev is quasi- 
homogeneous as well. Eqn. (76), or, equivalently, eqn. (77) can be used for recovering T* , and then for reconstructing 
S* by means of (18), if and only if one has ensured that SQ rev /f* is an exact differential, i.e., that SQ rev is integrable 
and /* is an integrating factor for 5Q rev . Otherwise, (76) is not even meaningful (a closed form on the left side of 
(77) should be equal to a non-closed form on the right side of the same equation). Cf. also [1]. 
Notice that, if S* has degree q, then (75) becomes 



i — 1 



,r-^ „ »_JL _ £ 2 | i_J_SL | . 0. (78) 



Q 

See also Appendix C. 



B. case where independent variables of weight zero appear 

Let us assume that, the independent variables X 1 * have weight zero for i = p + 1, . . . , n, with p < n. This means 
that £* have degree a for i = p + 1, . . . , n. The Euler vector field is 

D = aU*^-+Y a, X 1 *-^—, (79) 
dU* ^ dX % * v ' 

i<p 
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and the integrating factor is 

f* = aU*-J2 OigX"- 

i<p 

Then, one has 

i<p 

It is easy to show that the Gibbs-Duhem equation in this case is 

\ / i<p \ / p<i<n 

One can also obtain the analogous of eqn. (76) 

p<i<n 

(77) still holds. 

Concerning the Gibbs-Duhem equation in standard thermodynamics, see [47]. 



C. examples 

In the black hole case, as well known, one has [we don't write the asterisk in what follows] 

M $Q nj 



and 



S = 

2T 2T 



dS=^ dM -^dQ-^dJ; 



the generalized Gibbs-Duhem equation, associated with the quasi-homogeneity of S, is 

2T 2 \T ) 2T 2 \T \T 



which can be rewritten as [14] 



Then, one can find 1 /T from 



(f) j(- 2 dM+ l -Q 2 d(^)+j <m). 



It can be easily shown that 1/T can be determined apart from an undetermined multiplicative constant. 
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X. ENERGY REPRESENTATION AND LEGENDRE TRANSFORM 



In the energy representation, the fundamental equation is U* = U*(S*,X 1 *,..., X n *). One has 

dU* = T*dS* + C dX * * ■ ( 89 ) 

i 

U* is a quasi-homogeneous function of degree r and type (1, ct\, . . . , a n ); T* is quasi-homogeneous of degree r — 1, 
are quasi-homogeneous of degree at — r. (89) represents an exact quasi-homogeneous Pfaffian form of degree r. The 
Euler operator is 

D = S* t^+Y a t X l * -J— . (90) 
dS* ^ dX 1 * K ' 

i 

One has then 

U* = - (T*S* +J2 ai i* X 1 *). (91) 

i 

[Notice that the generalization to the case where S* is quasi-homogeneous of degree q is trivial] . 

As far as the Gibbs-Duhem equation is concerned, one easily finds [we shift to the case where S* has degree q] 

q s *-^+i d ^s^T*^ +J2~X l d(x 1 *^ £*) = 0. (92) 

For a proof see Appendix C. One can e.g. determine T* from 

d (s*^T*) = S* 3 ^- 1 -J2a t X* *~^f +1 dfx 1 *^ $ \ . (93) 

A formula analogous to (82) can be obtained if some independent zero-weight variables appear; one obtains (cf. 
Appendix C) 



S*~^r +1 d(s*Vr) +J2 — X 1 * "«' r+1 dfx'***''' $\ - &dX l * = 0. (94) 

i<p ^ ' p<i<n 



A. Legendre transforms 

Legendre transforming of the potentials involves the usual procedure. A noticeable difference, with respect to 
standard thermodynamics, is that, if, between the n + 1 independent variables, no variable of weight zero appears 
and no intensive dependent variable occurs, then one can iterate the Legendre transformations of the potential n + 1 
times and find a potential which is not identically vanishing. Moreover, the would-be intensive variables appear in 
the Euler operator involving the new independent variables. Let us e.g. consider the quasi-homogeneous free energy 

F*(T*, X 1 *,..., X n *) = U* -T* S*. (95) 

It is a quasi-homogeneous function of degree r with respect to the Euler operator 

^(r-^^ + E^*^- (96) 

i 

One can Legendre-transform also with respect to the remaining n variables X 1 *, . . . , X n *: 

n 

G r * l+1 (T*, fl, . . . , C) = U* - T* S* - xi * £* # 0. (97) 

1=1 
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The corresponding Euler operator is 



D n+1 = ( r - 1) T* A + £ ( r _ «i) C ^. (98) 

i * 

Notice that, all the potentials which are obtained by Legendre-transforming the internal energy U* are quasi- 
homogeneous of degree r with respect to the corresponding Euler operator; all the Massieu-Planck potentials, which 
are obtained by Legendre-transforming the entropy S*, are quasi- homogeneous of degree one. In other terms, the 
Legendre transform preserves the degree of quasi-homogeneity (but it changes variables and weights: the function 

h*(£*,Y*) =g* -C X*. (99) 

which is obtained by a Legendre transform with respect to a variable £ of weight a of a function g(X* , Y*) which is 
quasi-homogeneous of degree r and weights (1 — a, i — /3), has weight (a, r — /3)). For a proof, see Ref. [9] (theorem 
2, Appendix A therein). 

1. case where weight-zero independent variables appear 

Let us consider the case where a variable rf of degree r and a variable I* of weight zero appears: 

dU* = T*dS* + T dX* + rf dl*. (100) 

The Euler operator is 

One can Legendre-transform three times, obtaining 

G 3 {T*,C,V*) = \ [(1 - r)T*S* - (r - (3)CX* - rrj*!*} (102) 



and 



As = (r - 1) T* JL + {r - (i) C A + rrj*^. (103) 



2. case with would-be intensive variable of degree zero 



Let us consider the case where a variable (* of degree zero and a variable Z* of weight r appears: 

dU* = T*dS* + C dX* + C dZ*. (104) 

The Euler operator is 

D-S-JL + pX'^+rC-JL;. (105) 

One can Legendre-transform three times, obtaining 

G 3 (T*, r, C) = - [(1 - r)T*S* - (r - j3)^*X*] (106) 
r 

Ds = (r-1)T* ±- + {r-(3)e A (107) 

This case is analogous to what happens in standard thermodynamics e.g. when one passes from U(S,V, N) to 
F(T,V,N) and then to G(T,p,N). 



and 
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B. from the entropy to the energy representation and further considerations 



In order to pass from the entropy representation S*(U*, X 1 *,..., X n *) to the energy representation 
U*(S*, X 1 *,..., X n *), one inverts the first relation with respect to U* , which is possible because dS*/dU* > 0. 
If S* has degree one and weights (r, ct\, . . . , a n ), then U* has degree r and weights (1, a\, . . . , a n ). Analogous con- 
siderations hold when passing from the energy to the entropy representation. We show that, when it is possible to 
invert, at least locally, a quasi-homogeneous function with respect to one variable, a quasi-homogeneous function is 
obtained again, having obvious degree and weights. 

In the general case, let us consider a quasi-homogeneous function w = g(x x , . . . , x n ) of degree r and weights 
(ai, . . . , a n ); the partial derivatives 

are quasi- homogeneous functions of degree r — ai for i = 1, . . . ,n. Let us assume e.g. that pi ^ 0. Then, at least 
locally, one can invert w with respect to x 1 and find x 1 = h(w,x 2 , . . . ,x n ) such that g(h(w,x 2 , . . . ,x n )) — w. The 
inverse function h (where it exists and it is sufficiently smooth) is easily shown to be a quasi-homogeneous function 
of degree a\ and weights (r, 02, . . . , a„). In fact, one has 

g(\ ai x\ \ a2 x 2 , \ a ~x n ) = g(X ai h(w, x 2 ,..., x n ), \ a2 x 2 , . . . , A a "x") = X r g(x\ . . . , x n ). (109) 

Moreover, one has 

g(h(X r w, \ a2 x 2 , A a "x"), X a2 x 2 , . . . , A a "a; 11 ) = X r w = X r g(x\ . . . , x n ), (110) 

thus it has to hold 

h(X r w, X a2 x 2 , . . . , X a "x n ) = X ai h(w, x 2 , . . . , x n ). (Ill) 

Cf. Ref. [35], where an analogous property is shown in passing from the entropy representation to the energy repre- 
sentation in standard (homogeneous) thermodynamics. See also [48]. 

As a further example, let us consider Kerr-Newman black holes. By inverting (where possible) the relation 
T = T(M, Q, J) with respect to M, one obtains (at least formally) M = M(T, Q, J). We know that M and Q have 
degree 1/2, J has degree one and T has degree —1/2. In the new variables, M — M(T, Q, J) is a quasi-homogeneous 
function of degree 1/2 and weights (—1/2,1/2,1). The Euler operator corresponding to the independent variables 
(T, Q, J) is 

15 15 d 



and it holds 



DM = ^ M. (113) 



XI. FURTHER INSIGHTS FROM BLACK HOLES 



A feature of quasi-homogeneous thermodynamics, shared also with homogeneous thermodynamics, is that in the 
expression of the metrical entropy and of the absolute temperature a undetermined multiplicative constant appears 
(cf. also [1]). Such a constant is phenomenologically fixed once the absolute scale of temperature is chosen. From 
the point of view of the actual analytic expression of this constant, it should be furnished by a statistical mechanical 
calculation. The black hole case can be again useful. By dimensional analysis, if we make the hypothesis that the 
entropy of a black hole can depend only on ft, c, G, kb, we find that 

S^^A; (114) 
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this is the same procedure as in the original proposal of Bekenstein [49] . For the temperature one obtains 



he ( OA \ 1 . . 

Toc ~h \TflGMWi) ^ < 115) 

A dimcnsionless constant has still to be determined. 

If one considers standard thermodynamics cases, this hypothesis which allows to recover the dimensions of the entropy 
by involving only the above constants works for the photon gas but it does not work in other cases; e.g. in the case of 
a Fermi gas of electrons, the electron mass represents a further scale, related to a microscopic analysis of the system, 
to be taken into account. In our case, the confirmation of the above hypothesis comes from the Hawking effect, which 
gives us the behavior of the black hole temperature, to be compared with (115). Then, because of the absence of a 
particle mass scale and, e.g. in the Schwarzschild and in the Kerr cases being involved only vacuum solutions of the 
Einstein equations, one could be tempted to think about some sort of "graviton gas", but, actually, this is far from 
being evident. See also Ref. [14]. 

Black holes represent an example of the the procedure of reducing/enlarging consistently the thermodynamics space 
discussed in sect. III. In fact, the construction of the fundamental equation sketched in sect. IV is allowed also in the 
case of a black hole with J = and of a black hole with Q = 0. Two parameters appear, M, Q and M, J respectively. 
The integrability condition is then trivially satisfied, but a well-defined potential S is obtained everywhere. Moreover, 
black holes of the Kerr-Newman family are a consistent extension of the black holes of the Rcissner-Nordstrom family 
(by adding J) or of the Kerr family (by adding Q). To this family belong also a thcrmodynamically degenerate case: 
the Schwarzschild case J = = Q, which is described only by one variable M. Notice also that the Kerr-Newman 
family can be obtained by putting A = in the KN-AdS family. 

In the case of black hole thermodynamics, it is tempting to conjecture that the quasi-homogeneous behavior of the 
thermodynamic potentials could be related to an explanation, to some extent, analogous to the one for the quasi- 
homogeneous behavior of the thermodynamic potentials in standard thermodynamics near the critical point [50,51]: 
As in standard thermodynamics the (leading order) quasi-homogeneous behavior can be related to the conformal 
invariance of the underlying quantum field theory near criticality [51] , one could conjecture that the quasi-homogeneous 
behavior of black hole thermodynamics could be related to some sort of conformal invariance of the quantum theory 
underlying General Relativity (superstring theory?). This field deserves further investigation. We however point out 
that it does not seem to be necessary that quasi-homogeneous thermodynamics should be related to some sort of 
conformal invariance in every case. For black holes, as well as for the general case, one should determine the reason 
of a quasi-homogeneous behavior by means of statistical mechanical calculations. In the following section, a further 
discussion is found. 



XII. HEURISTICS: QUASI-HOMOGENEOUS THERMODYNAMICS AND THERMODYNAMIC LIMIT 

In this section we suggest that the thermodynamic limit, under suitable hypotheses, can lead only to quasi- 
homogeneous thermodynamics. The thermodynamic limit should not be intended in a literal sense, of course, but as 
a tool which allows to determine the leading terms in the statistical mechanical functionals as large systems occur. 
Cf. e.g. Ref. [52], chap. 4, and also Ref. [53]. Thermodynamics of macro-systems emerges as an asymptotic law 
which is extrapolated from the asymptotic behavior of statistical mechanics. To some extent, the thermodynamic 
limit is also a tool by means of which the vague notion of "macroscopic system" is meant to be implemented, and 
a proper thermodynamic behavior recovered. Such a behavior depends on the interactions and on the nature of the 
system involved (e.g., fcrmionic or bosonic matter), together with a dependence on the initial conditions [from a 
quantum- mechanical point of view, a sort of "preparation" of the system, even if it is of "astrophysical size"]. As 
far as the gravitational interaction can be neglected, the standard thermodynamic behavior emerges, with the usual 
asymptotic laws which justify the extensivity of standard thermodynamics. At the same time, extensivity cannot be 
considered as an absolute and unmodifiable property of thermodynamics, because of the purely attractive nature of 
gravity. One can also expect a different asymptotic behavior for different systems at different scales. See also the 
discussion in the following subsections. 
The basic set (B) of assumptions is the following: 

Assumptions (B): 

Let {Zi} be the set of the thermodynamic variables, both independent (Z^p) and dependent (z\ D ^) ones. Let the 
independent variables describe the Gibbs thermodynamic space. 

Let A be the parameter such that, by taking A — ► oo, one recovers the thermodynamic limit. 
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Assume that the thermodynamic limit properly describes the bulk properties of the system under consideration. 
Assume that the thermodynamic limit of the statistical mechanical quantities, like e.g. the entropy or the free energy, 

exists under the hypothesis that, in the limit A — > oo, the independent thermodynamic variables z\ I \\) are rescaled 
as 

zf)=^(A)zf)°, (116) 

where <?j-^(A) is a positive function which is either an invertible function of A or it is one (or a constant), and 
is the value of the i-esime variable Z^p at an arbitrary reference state (cf. Rcf. [53] for the extensive case). 

We now state the following result (R-group). 

Result (R-group): 

Let us assume the set of assumptions (B) and that for the rescalings of the independent variables the group property 
holds: 

9?\\»)=g?\\)g\ I \»). (117) 

This implies that for each independent variable there exists a real number such that g^\X) = A Q * ' . 
Let us assume that the dependent thermodynamic variables 

Z[ D) =Z\ D \Z[ I \...,ZP), (118) 

are obtained as leading terms in the asymptotic expansion of the corresponding statistical mechanical functionals 
z\ D ^ in the following sense: for each i there exists a positive continuous function Pi(X) such that 

lim *W(A) Z[>\ . . . , 5 (/>(A) ZW) = . . . , Zi% (119) 

A^oo Pi(A) 

and that the domain of the functionals z\ D ^ is invariant under the above rescalings. 
Then, the dependent thermodynamic variables are quasi-homogeneous, i.e. they satisfy 

Z\ D \£\\) Z[>\ . . . , 5 W(A) ZW) = 9 l D \ X ) ZP\Z['\ ZW) (120) 
where g\ D] '(A) = A Q ^' and a| D) G IR. 

Moreover, for all i, the functions p,(A) are regularly varying functions [54]: 
for all p e IR+ 

lim ^=^ D '. (121) 

A^oo Pi(X) 

(a| D ' is also called degree of pi) . 

In order to be more explicit, let us consider (R-group) for a specific case. Let us e.g. assume that the thermo- 
dynamic limit for the statistical mechanical functional S representing the entropy exists in the sense that 

lim — ^— S(\ au U, \ av V, \ aN N) = S(U, V, N), (122) 

A^oo ps{\) 

where the function ps is positive and ajj, ay, o>n are real numbers. Then for each real p > the function ps satisfies 

Um ps(Xp)_ = ^ (123) 

A^oo ps(\) 

for some 7 G IR. Moreover, S is quasi-homogeneous of degree 7. In fact, from (122) it follows 

lim — tz — r <S((A p) au U, (A p) av V, (A p) aN N) = S(U, V, N). (124) 

A^oc p s (\ jj.) 
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Moreover, it holds 



lim — — S(X au (p au U),X av {p av V),X aN (p aN N)) = S(p au U,p av V, p aN N). (125) 

Then 

ps{\) S{{\ p) a " U, (A p) a - V, (A p)<*» N) 



A^+oo p s (X p) 5(A Q " (^u u),X a v (fi a ^ V),A Q « (p aN iV)) 

5 ^' V '^ - Km ->#> (126) 



S(p au U,p av V,p aN N) x^oo p s {Xp)' 

which is possible only if the above conditions are implemented. A general proof is found in Appendix E. Cf. [55], 
where a rigorous approach for asymptotics of tempered distributions is developed. 

The case where one or more variables are "intensive", i.e. are not rescaled, presents no difficulty The above rea- 
soning about the leading order in the asymptotic law identified in the thermodynamic limit is corroborated by the 
standard thermodynamic limit, which allows to find a leading-order homogeneous thermodynamics. We define quasi- 
homogeneous thermodynamic limit the thermodynamic limit satisfying the group property. 
Moreover, we notice that, if e.g. for the thermodynamic entropy S(U, V, N) it holds 

S( 9u (X) U Q , 9V (X) V ,g N (X) N ) = g s (X) S = g s (X) S(U ,V ,N ), (127) 

where the functions gu, gv, 9n, gs are positive and the group property holds, then S has to be quasi- homogeneous 
and the are power-like functions. The key point is that a behavior like the one in (127) under the group property 
is necessarily quasi-homogeneous; as it is shown in Ref. [9], "powers need no generalization" [see Appendix A in Ref. 
[9] and Appendix C in Ref. [10]]. A different proof is sketched in Appendix D for the sake of completeness. The 
aforementioned theorem implies that actually one has 

S(g(X) a " U,g(X) av V, g(X) aN N) — g(X) as S(U,V,N) (128) 

with the same scaling function g(X). By defining 

A = 5 (A) (129) 

one finds S(X au U, X av V, X aN N) — X as S(U, V, N), i.e. the standard definition of quasi-homogeneous behavior for a 
function. A generalization of the quasi-homogeneous behavior of a function under rescaling with the group property 
is not allowed. 

We now relax the requirement for the group property, one can still make a conjecture (C-reg) which would allows 
a quasi-homogeneous behavior under rather general conditions. 

Conjecture (C-reg): 

Let us assume the set (B) and that g, gi (for i = 1, . . . ,n) are positive and continuous functions which are regularly- 
varying functions of real degree 7, ai, . . . , ct n respectively. 
Assume that the statistical mechanical functional Z satisfies 

lim -}-Z(g 1 (X)Z 1 ,...,g n (X)Z n ) = Z(Z 1 ,...,Z n ), (130) 
where Z l are the independent thermodynamic variables (we omit the suffix (I)). 

We conjecture that, under suitable hypotheses, the asymptotic Z{Z X , . . . , Z n ) is quasi-homogeneous of degree 7 and 
weights (ai,.. . ,a n ). 

For each positive real a< one should find 

lim —L-Z( gi (Xp)Z\...,g n (Xp)Z n )=Z(Z\...,Z n ) (131) 

A^oo g[A p) 

" = " lim _ L-Z(. 9l (A) M ai Z\...,g n (X)p a " Z n ) (132) 

A^oo g{\) pi 

= ± Z(p a - Z\...,p a " Z n ) (133) 
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where "=" indicates that further hypotheses on the statistical mechanical functional should be satisfied. Notice that 
the case gi(X) = A Qi is a sub-case of the preceding one and requires only that the limit (130) exists and that the state 
(fi ai Z 1 ,. .. , yu"" Z n ) belongs to the domain of the statistical mechanical functional Z. 

From a physical point of view, it can be noticed that, when passing from the statistical mechanical functional S 
representing the entropy to its asymptotic S sm under dilatations, as in the standard thermodynamic limit, it would 
be not strictly necessary, a priori, to impose that g(X) — X in order to obtain an homogeneous S sm [herein we indicate 
by S sm the asymptotic which is indicated by S in (122). S sm corresponds to the statistical mechanical entropy one 
takes into account in sect. IIIC2]. A generic regularly varying function of degree one could also be allowed (e.g., 
g(X) = X log(A) could be allowed). The requirement 

lim i S(X x) = S sm (x) (134) 
amounts to the requirement that in some sense also S is extensive. 

We have found that, under rather general hypotheses, the only possible outcome of the thermodynamic limit is a form 
of quasi-homogeneous thermodynamics. The thermodynamic limit, at the same time, is not expected to describe all 
possible systems; e.g., it cannot be applied to a thin (molecular thickness) metallic film [56]. Nevertheless, it can be 
safely applied to a huge class of macroscopic systems, and, in the non-relativistic case, it allows to conclude that the 
Thomas- Fermi approximation becomes exact in the limit. 

In partially concluding these heuristic considerations, we point out that a relation between a so-called pseudo-extensive 
thermodynamics and the thermodynamic limit has been postulated in Ref. [57] and an analysis of Newtonian self- 
gravitating systems is made in Ref. [58]. We limit ourselves to refer the reader to the aforementioned papers. 

A. non-conventional thermodynamic limit 

A non-conventional thermodynamic limit has already been proposed. A very interesting discussion of the thermo- 
dynamic limit in physical systems is found in Ref. [59] . Therein, physical systems are divided into two classes each of 
which can be associated with its own thermodynamic limit; one can distinguish between the standard thermodynamic 
limit (STL) and the inhomogeneous mean field thermodynamic limit (IMFTL) which occurs e.g. in self-gravitating 
systems as the one of HNT model. The latter kind of thermodynamic limit can be related with the presence of 
long-range forces [59]. In order to identify the correct limit, the key notion of characteristic bulk length scale is 
introduced, which corresponds to a typical length scale X typ resulting from the interaction of many particles in the 
system. Such length scale Xt yp plays the role of characteristic invariant in the thermodynamic limit. In the case of a 
classical self-gravitating isothermal gas, the Jeans length Xj ean s is such a typical length. Moreover, the requirement 
that mean thermodynamic quantities exist almost everywhere in the thermodynamic limit is a tool for selecting the 
appropriate thermodynamic limit at constant Xt yp [59]. A complementarity between STL and IMFTL can also be 
allowed for the same system, in the sense that STL and IMFTL can describe complementary asymptotic properties for 
the same large but finite system (the case of classical Coulomb systems is discussed) . We don't discuss further on this 
approach herein, we limit ourselves to point out that one obtains the same limit, in the case of classical self-gravitating 
matter, if one keeps Xj ean s fixed or £ fixed (cf. sect. V). Notice that £ is an intensive function with respect to the 
Euler operator relative to the HNT model. At an heuristic level one expects that in general there exists a physically 
meaningful function which is left invariant under a quasi-homogeneous thermodynamic limit and which is intensive 
under the corresponding Euler operator. 

A non-conventional discussion of the thermodynamic limit is found also in sections 1.2 and 4.2 of Ref. [27]. Both for 
stable and for unstable interactions are discussed, and the criterion for selecting the limit is a comparison between the 
kinetic and the potential energies involved in the system. According to [28] [p. 1-3] , a non-conventional thermodynamic 
limit has to be associated with the peculiar scaling behavior for "non-stable interactions" ; the notion of "Thomas- 
Fermi thermodynamic limit" is also introduced for non-extensive systems. Note that a more general thermodynamic 
limit can also be allowed, where also the coupling constants are allowed to vary in the limit. Cf. [25,28,59]. 

B. absolute weights again 

It is evident that one has a reference quantity with respect to which all the others are "weighted" in the thermody- 
namic limit. Then an overall undetermined scale factor q is expected to appear, unless an absolute weight is known. 
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In the standard homogeneous case, the problem of assigning absolute weights does not arise, they are all equal to one 
for extensive variables and zero for intensive ones, the thermodynamic limit is performed with the aim of studying bulk 
properties of the system, and the ratios U/N and V/N are kept constant. Roughly, one fixes the energy per particle 
and the volume per particle (or, equivalently, the density N/V and the energy density U/V) and then performs the 
limit with the aim of neglecting any finite-size effect. The same is true also for a massless particle gas, like the photon 
gas, where the limit U, V — ► oo, U/V = const, is performed. 

If a strict quasi-homogeneous picture is required, then the above problem is instead present. When N is available, as it 
is when non-zero rest-mass particles are considered, it is to some extent natural to assume that its weight is one. Such 
a special role of iV can be justified. In the model of Hertel, Narnhofer and Thirring, the existence of the Boltzmann's 
entropy per particle in the thermodynamic limit liniTv^oo log(f2)/AT is ensured, in agreement with Boltzmann's pos- 
tulate, which states that the entropy per particle \og(Q)/N exists for an equilibrium system. As a consequence of the 
coincidence of liniTv^oo log(f2)/JV with the thermodynamic entropy density (per particle), one can justify a posteriori 
the assumption that N has weight one, in fact log(fi) in the thermodynamic limit coincides with the thermodynamic 
entropy which has degree one. The other weights a» appropriate to the physical condition one is considering for 
the system are recovered in the asymptotic limit. Then, no ambiguity for an overall multiplicative constant in the 
definition of the weight like the one discussed in sect. IIIC appears. Note also that the Boltzmann's entropy being 
the logarithm of the number of micro-states compatible with a given macro-state, by involving the (vague) concept of 
macro-state and thus of macroscopic system can legitimately related to the thermodynamic limit (again in the sense 
of asymptotic law). [About Boltzmann's entropy and its coincidence with the Gibbs entropy (which, in the quantum 
case, corresponds to the von Neumann entropy) see e.g. [60].] A more puzzling case is represented by a self-gravitating 
massless gas, but weights can be still fixed. See sect. VI. 

In general, if A — > oo corresponds to the thermodynamic limit, one could also require the existence of the Boltzmann 
entropy density functional as A — > oo, to be defined as 

s bo it z = lim — ^-log(fi), (135) 

A— >oo A s 

where as is a constant which is fixed by the requirement of finiteness of the above limit, as has to be introduced 
because the scale A can also be such that A/A is not asymptotically a positive constant, i.e., A ^ N in the sense of 
the asymptotic behavior. This more general choice for the scale is also necessary if N is not a good thermodynamic 
variable (as in the massless gas case). We know that it is consistent to fix the weights by choosing the Boltzmann 
entropy to be of degree one. 

If another variable, e.g. the internal energy, would be used as reference variable, the relative weights could be still 
recovered in line of principle, but the aforementioned ambiguity would emerge in absence of a criterion for finding the 
absolute weight at least of one (independent or dependent) variable. 

One could also choose to work in the energy representation, and decide to fix the weight of the internal energy to 
be one. Then the entropy in general would not be a degree one quasi-homogeneous function; in the HNT model, one 
would obtain deg(V) = —3/7; dcg(TV) = 3/7; deg(5) = 3/7, against the weight 7/3 of U resulting by assuming that 
N has weight one. This procedure, although legitimate because of the overall ambiguity, is less satisfactory, because 
the thermodynamic construction, in the case the entropy is a quasi-homogeneous function and the weights present 
the overall factor ambiguity, allows naturally the thermodynamic entropy to be of degree one. The only exception 
which, in line of principle could be admitted, is given by the case where the weights in the Gibbs space are all fixed 
and the statistical-mechanical entropy is not of degree one, but it is apparently not allowed (unless some reason for 
fixing the absolute weight of a variable appears). 

C. thermodynamics, weights and gravity 

As far as gravity is negligible, one finds the standard extensivity property for the fundamental equation in the 
entropy [energy] representation; when the scale of the system is such that is no more possible to neglect self-gravity 
effects, then the weights of the thermodynamic variables change according to the nature of the system and to the scale 
itself. The quasi-homogeneous character of thermodynamics with fixed weights (a\, ...,a n ) has to be intended as 
relative to the validity range of the asymptotic expansion which leads to those particular (ai, . . . , a n ), whose nature is 
thus non absolute. [One could also be tempted to write = aj(A), in order to recall that the weights can change with 
the scale. Nevertheless, in light of the above theorem about quasi-homogeneity, such a dependence is not allowed. The 
quasi-homogeneity is extrapolated at each fixed scale, and it is made "absolute" by the thermodynamic formalism.] 
Also an estimate of the error for finite A should be allowed [52]. One could think also to a sort of "evolving picture" 
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for the formalism, in which one starts from standard homogeneous systems, characterized by stable interactions, and 
then, as far as gravity becomes non-negligible, a quasi-homogeneous behavior is allowed. Weights can change, even in 
a discontinuous way, depending on the physical conditions leading to the new stationary equilibrium state. Because 
of the Hawking effect, the "end of the story" does not occur even when the scale and the physical conditions are such 
that a black hole forms, whose "classical" weights are the ones presented in sect. IV. 

It is interesting to underline again that ordinary thermodynamics implements the picture described above, as discussed 
in several fundamental papers on the stability of matter [61,62,29,63]. See also Refs. [64,43]. Strict quasi-homogeneity 
is expected to be effective when gravitation begins playing an essential role in the physics of the system. Moreover, 
the so-called fourth law of thermodynamics [13] should be generalized in such a way to include the general quasi- 
homogeneous behavior of the thermodynamic variables [dependent as well as independent ones] . 

XIII. CONCLUSIONS 

Quasi-homogeneous thermodynamics is proposed as the unifying picture for thermodynamics where both stationary 
black holes, fcrmionic non-relativistic matter, and self-gravitating electromagnetic radiation are found. 
We have chosen to introduce quasi-homogeneity in thermodynamics by generalizing the formalism developed for 
standard homogeneous thermodynamics [1]. The basic requirement is that the Euler vector field which generates 
quasi-homogeneous dilatations in the Gibbs' thermodynamic space is a non-trivial symmetry for the integrablc Pfaffian 
form 8Q rev . This means that the symmetry maps leaves of the thermodynamic foliation onto other leaves. The leaves 
of the foliation are naturally associated with a quasi-homogeneous function S* of degree one, which is suggested to 
be the (unique apart from a multiplicative constant) metrical entropy corresponding to the fundamental equation 
in the entropy representation. This is true also in the case one knows the weights of the thermodynamics variables 
apart from an overall unknown factor. Notice that quasi-homogeneity could be introduced also at the level of the 
fundamental equation in the entropy representation or in the energy one, but, in the case where only the relative 
weights of the independent thermodynamic variables are known, it should be considered equivalent the choice of a 
degree one internal energy with respect to the choice of a degree one entropy (cf. also [9]), whereas the aforementioned 
approach privileges a degree one entropy. 

The consequences of the lack of homogeneity arc recalled, and superadditivity of the metrical entropy is a privileged 
property ensuring the second law even in absence of concavity for S* . 

The generalized Gibbs-Duhem identities are also derived, and the energy representation and the behavior under 
Legendre transforms is analyzed. 

An heuristic argument which can relate quasi-homogeneous thermodynamics to the thermodynamic limit in statis- 
tical mechanics (even if in a non-conventional framework) has been introduced. According to this argument, under 
the group property only a form of quasi-homogeneous thermodynamics can be recovered in the thermodynamic limit. 
Moreover, a mean field thermodynamic limit appears to be appropriate; the suggestion from the existing models is 
that one has to search for asymptotic scaling properties of the equations allowing to find out the mean field ther- 
modynamics for the system, in fact the scaling properties of a finite temperature Thomas-Fermi equation appears to 
be involved in the quasi-homogeneity of the HNT model, whereas the scaling properties of the TOV equations are 
related with the quasi-homogeneity of self-gravitating radiation. 

Gravity appears to play a fundamental role in allowing a generalization of the standard scheme for thermodynamics. 
The purely attractive nature of gravity, as well-known, is the reason for the failure of the extensivity property in 
thermodynamics of self-gravitating systems, because of the absence of saturation which leads to an implosive insta- 
bility, and matter systems eventually can implode into a black hole if the mass is sufficiently large. Moreover, the 
thermodynamic ensembles become inequivalent [28], and the lack of concavity leads to phase transitions in matter 
systems in the canonical ensemble [27,28]. Negative heat capacities arise in the micro-canonical ensemble. Also black 
holes present this feature. 

Black hole thermodynamics requires some special comments. From the point of view of thermodynamic formalism, 
black hole thermodynamics can be realized to be a rather straightforward generalization of the formalism developed 
for standard thermodynamics in [1]; moreover, its quasi- homogeneous behavior can no more be considered as excep- 
tional, because self-gravitating systems exist where quasi-homogeneous scaling laws are satisfied. Nevertheless, from 
other points of view black hole thermodynamics is still to be considered special, also because it is still unclear which 
statistical mechanics should lie beyond it. 

We limit ourselves herein to this proposal for a generalization of the thermodynamic formalism, without pursuing 
it further on [apart from a note on the third law in Appendix B]. For future investigations, it would be interesting to 
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develop a statistical mechanical formalism where General Relativity is self-consistently included. 
It would be also interesting to explore, if other long range forces are allowed to be included in the framework of 
quasi-homogeneous thermodynamics. Classical two-component Coulomb matter seems to belong to this framework 
[59] , further investigations are required in order to see if this is true also for quantum Coulomb matter. 
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APPENDIX A: INTEGRAL OF co/fi AND POTENTIALS ASSOCIATED WITH 
QUASI-HOMOGENEOUS EXACT PFAFFIAN FORMS 



We first present some results in the general case of a generic symmetry X for a PfafHan form u>. Then we translate 
the results for the case of the quasi-homogeneity symmetry. We recall that a vector field X is defined to be a symmetry 
for a Pfaffian form lo if 

(L x lu) Alu = 0. (Al) 
This means that there exists a function h such that 

L X Lo = huj. (A2) 

See [5,6]. 

We need the following two lemmas: 
lemma 1 

Let o) be an exact Pfaffian form; 
let X be a symmetry such that 

L x w = (A3) 
where q ^ is a constant. Then, one finds that uj = d<p is implemented by 

(f>=- i x u>. (A4) 

proof: The proof is elementary. One has trivially 

dxj> = - d(ix u>) = — ix(duj) H — L x u>- (A5) 

From dui — because u is closed and from the symmetry of ui, which implies L x to = q u> the thesis follows. Notice 
that 4> satisfies X<p = q <j). □ 



lemma 2 

Let W(o) be an exact Pfaffian form; 

let X be a symmetry such that L x lu^ = 0; 

let W be the associated potential: dW = W(o) . 
Then 

dF 

dW=— (A6) 

where F satisfies 

X F = q F, (A7) 

with q G IR constant. Moreover, 

ix W(o) = ?■ (A8) 
As a consequence, q — is allowed if and only if the symmetry is tangent (i.e. trivial). 

proof: We have 

= L x uj (q) = L x dW = d L x W. (A9) 
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As a consequence, 



L x W = X W = q (A10) 
where q is a constant. We can define a positive definite function F such that 

W = \og(F). (All) 
Then we get that (A10) is equivalent to the following equation for F: 

X F — q F; (A12) 



moreover 

dF 

dW = — , (A13) 

and 

ix W( ) = ix dW = L x W = 1 X F = q. (A14) 

If the symmetry is tangent, then ix W(o) = 0. Then g = necessarily. If q = 0, then and W satisfy If = and 
A = respectively. As a consequence, ^(o) = ix dW = Lx W = 0, which completes the proof. □ 
Notice that W has to be found by quadratures, the contraction of the Pfaffian form with the vector field X is not 
useful in order to find a potential without explicit integration. 

One can in general show that the following theorem holds: 
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theorem 1 



Let to be a integrable PfafEan form which is defined in a connected, simply connected domain; 
let X be a non-trivial symmetry for ui; 
let fi = ixu. 

Then n. is an integrating factor for u> and the foliation corresponding to the PfafEan equation u — is associated with 
a potential F which satisfies the following equation: 

X F = F. (A15) 

proof: The proof that oj/u. is closed is trivial. See also [6] and [1]. Moreover 

Lx - = 0, (A16) 

being X a symmetry for u. In fact, 

Lx — = 4j ((L x cj)fj, - (L x h)oj) = ,. 1 ,„ ({L x u>)(ixu) - (L x ixu)Ld) (A17) 

/i fi z [ixuy 

i x (luA(L x uj)), (A18) 



(ixu) 2 

where the rules ixLx = Lxix and ix(ce A (3) — (ixct) A f3 + (— )™ a A {ix/3) for a n-form a and a m-form f3 are used. 
Cf. [65]. 

As a consequence, uj/fj, is an exact one- form such that the previous lemma applies. One has 

-=dW=™ (A19) 

where F satisfied eqn. (A7) with q ^ 0, because of the hypothesis of non-trivial symmetry. Then one finds 

dF = — u, (A20) 
M 

which satisfies Lx dF = q dF and is closed. Then, from (A4), we obtain 

F =- — ixu;, (A21) 
q ii 

which, because of the definition of fx, implies q = l. □ 
Notice that the integrating factor satisfies 

Lx M = Lx ixu = ix Lx lu = h fi. (A22) 

It is interesting to underline that this theorem allows to construct a general theory of thermodynamics if 5Q rev has 
a non-trivial symmetry. It is easy to see that the constructions of sects. Ill A and III B remain unaltered, and also 
the constructive assumptions of sect. VIII can be trivially generalized [it is sufficient to delete "quasi-homogeneous" 
in al) and a2) and to change a5) into "the metrical entropy S satisfies X S = S]. The case of quasi-homogeneous 
symmetry is immediately obtained by requiring that each independent variable in the thermodynamic space has the 
same kind of symmetry as the potential S. More general symmetries can also be introduced. Further details will be 
given elsewhere [66]. 

1. the quasi- homogeneous case 

The following results are immediate consequences of the previous lemmas and theorem in the case of quasi- 
homogeneous symmetry. 

result 1: 
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Let 



n 

UJ = 

i=l 



Ui dx l (A23) 



be an exact quasi-homogeneous Pfaffian form of degree r ^ and weights (ai, . . . , a„); 
let 

n ft 

Ds T, a *^w (A24) 

i=l 

be the Euler operator; 

let the symmetry be nontrivial (or transversal), i.e. 

n 

u>(D) = i D uj = ^ a i uj i x i ^0; (A25) 
Then, one Ends that to — d<p is implemented by 

<f>=- uj(D). (A26) 
r 

The above result can be used even when some zero weights occur. It cannot be used if r = 0, in which case the 
potential can be found only by quadratures. When one considers e.g. SQ rev /f, the potential has to be found by 
integrating the exact form. 

result 2 : 

Let W(o) be an exact quasi-homogeneous PfafEan form of degree r = 0; 

let W be the associated potential: dW = W( ) . 
Then 

dF 

dW=— (A27) 

where F is a quasi-homogeneous function of degree q with respect to the Euler vector field D. Moreover, 

iD W( ) = q. (A28) 
As a consequence, q — is allowed if and only if the symmetry is tangent. 

Notice that W has to be found by quadratures, the contraction of the Pfaffian form with the Euler vector field 
is not useful in order to find a potential without explicit integration. 

result 3 : 

If u>f r ) is a quasi-homogeneous integrable Pfaffian form of degree r defined in a connected, simply connected do- 
main, if the symmetry is non-trivial and if [i is the corresponding integrating factor to(D), then 

d W^^l = d 4 (A29) 
jj, F 

where F is a quasi-homogeneous function of degree one with respect to the Euler vector field D. 

The degree r of W( r ) can also be zero in the case where at least one weight, say ai, is different from zero (the 
would-be intensive variable u\ has then weight — ai). Herein examples are given where < r < oo. It would be 
interesting to see if r < is physically meaningful. 

Notice that the fact that the case q = is not allowed by the requirement of non-triviality for the symmetry can be 
verified by direct inspection, in fact from eqn. (A20) it follows 
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d i F= — w i Vi = l,...,n, (A30) 
M 

thus DF = becomes 

" F 

^ a t x l = 0, (A31) 

which is not allowed, in fact, because of [i = Y^i=i ai x * w *' ^ would imply F = 0. It is interesting to note that, from 
a mathematical point of view, W is an almost quasi-homogeneous function of degree zero with first order deficiency 
index identically equal to one (cf. [3]), i.e. its behavior under quasi- homogeneous scaling is 

W(X ai x\ X a "x n ) = Wix 1 , ...,x n )+ log(A). (A32) 

This property of W is to be related with the non-triviality of the symmetry, which allows to find an integrating factor 



2. some more results on quasi- homogeneous integrable Pfaffian forms 

We have found that, for non-trivial symmetry, 

W( r ) = fif(r-i) dF, (A33) 
where 5( r -i) is quasi-homogeneous of degree r — 1. If G is a quasi- homogeneous function of degree q such that 

\r) = 9(r-q) dG (A34) 



then necessarily 

where Q — const. In fact, one has 
and 



dG/dF 

Moreover, the quasi-homogencity of G(F) implies 



G = C F q , (A35) 
W( r ) = 9{ r -q) dG = g {r _ 1} dF (A36) 



9ir- q) = (A37) 



An 

D G = q G = — F, (A38) 



dF 
that is, 

dG G 
dF=«F> 

whose solution is (A35). For q = 1, one obtains that the quasi-homogeneous function _F of degree one implementing 
(A33) is unique apart from a multiplicative constant. [Notice that this result can be generalized easily to the case of 
a generic symmetry X. By referring to theorem 1, if G satisfies X G = q G and u = b dG, then G = G(F), where F 
is the potential (A15), and G — k F q , where n is a constant]. 

It is evident that, if D is a symmetry generator, D q = q D is another symmetry generator for q G IR — {0}. It is 
also interesting that, by changing q, the integrating factor f q changes in such a way that uj/f q = dlog(F q ), where 
F q = F 1 ^ is quasi-homogeneous of degree one with respect to D q (i.e., D q F q = F q ). In fact, one has f q = q f and 

f = lj = l dl °S( F ) - dlogiF 1 ^); (A40) 



then, D q F x li = qD F x l q =ql/q F 1 ^' 1 D F — F 1 /*. 
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APPENDIX B: S AS A FUNCTION OF T IN QUASI-HOMOGENEOUS THERMODYNAMICS AND 

THE THIRD LAW IN ENTROPIC FORM 

We assume that the metrical entropy S* is quasi-homogeneous of degree one, and that a, a±, . . . , a n are the weights 
of the independent variables U* ,X X * , . . . ,X n *. The temperature is then a quasi- homogeneous function of degree 
a— I. Let us consider S* — S*(T*, X 1 *,..., X n *). The entropic form of the third law, for standard thermodynamics, 
states that S —* So as T — > 0, where Sq is a constant which has to be zero [1]. Planck's restatement of the third law 
S — > as T — > is mandatory, if the third law holds, also in the case of quasi-homogeneous thermodynamics. In fact, 
if, for any fixed value of X 1 *,..., X n *) 

Hm S 1 * (T*, X 1 *,..., X"*) = S*{0,X 1 *,..., X n *) = S* (Bl) 

holds, then, from 

lim S*(X a - 1 T*, X^X 1 X a "X n *) = X lim S*(T*, X 1 *,..., X n *) = X SX (B2) 

and 

lim S*(\°- 1 T*, X^X 1 *,..., X a "X n *) = S*(0, X^X 1 *,..., X a "X n *) = S* (B3) 

one finds Sq = 0. In the above formulas, continuity of S at T = is assumed (otherwise the third law is violated 
[67]). 

We know that black hole thermodynamics violates the entropic form of the third law. This behavior is not a general 
feature of quasi-homogeneous thermodynamics. We construct a toy-entropy which is both quasi-homogeneous and 
superadditivc and implements the third law. Herein, as in sect. VII, with x we indicate collectively the n variables 
appearing in S. Our starting point consists in realizing that, given two non- negative superadditive functions g(x), h(x), 
the function 

F(x) = g(x) h(x) (B4) 

is superadditive too. In fact, 

F(x! + x 2 ) = g(xi + x 2 )h(x 1 + x 2 ) > (g(xi) + g(x 2 ))(h(x 1 ) + h(x 2 )) (B5) 
= g{x 1 )h{xi) + g(x 2 )h(x 2 ) + g(x 1 )h(x 2 ) + g(x 2 )h(xi) (B6) 
> gfaMx!) + g{x 2 )h{x 2 ) = F( Xl ) + F(x 2 ). (B7) 

Then, let us consider g = 2y/U*X* and h = X* 2 . Both are superadditive functions (g is superadditive because of 
its homogeneity and concavity). Then S* = 2\/U*X* X* 2 is superadditive and quasi-homogeneous. For S* as a 
function of T* , one finds 

S* = 2 T*X* 5 (B8) 

which satisfies S* — > as T* — > 0. 

It is also interesting to notice that, if S is a strictly superadditive function everywhere (i.e., also on the boundary 
T = 0) then the entropic version of the third law in Planck's form has to be violated; in fact, if (Uo, Xq, . . . , X$) is 
such that 

T(U ,X \...,XZ) = 0, then 

S(U + U ,X 1 +Xt,...,X n + XS) > S(U a , X \..., X%) + S(U, X u ..., X n ) (B9) 

is possible with the strict inequality only if S(Uq, Xq, . . . , Xq) > 0. In the case of a self-gravitating system belonging 
to quasi-homogeneous thermodynamics framework, if the energy is strictly subadditive everywhere, then the third 
law cannot hold. 

As far as mathematical properties of the Pfaffian form SQ rev ensuring the validity of the entropic form of the third 
law are concerned, it can be shown that a quasi-homogeneous 5Q rev which is C 1 everywhere is a sufficient condition, 
in analogy with standard thermodynamics [1,67]. In fact, a superadditive non-negative function S : T> — > IR + U {oo} 
cannot be divergent at a limit point xq of its domain [we mean that all the points x in the domain satisfy ||x|| < oo, 
where || • || stays for the Euclidean norm and that xq can be included in the domain by putting S(xq) = oo]. This is 
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evident if one translates rigidly the domain in such a way that xq = 0. Then, if zq is another point in the domain of 
S one has S(z ) — S(z n + 0) > S(0) + S(zo) = oo. Thus, if S ^ oo, then S has to be finite. Then S has to be finite 
in the limit as T — > 0. As a consequence, if / is the integrating factor, one has 

df dT* 
' =1 + 5*^ (BIO) 



dU* dU* 

and dT* /dU* = 1/C, where C is the heat capacity at fixed parameters X % *, which has to tend to zero if S* is finite 
and non- negative as T* — > 0. Then, if / is everywhere C 1 (and it is such if SQ rev is C 1 everywhere), it holds S* — ► 
as T* — > 0. Moreover, it is evident from (18) that lim T »^o S* = holds if and only if J 5Q rev /f* —> -oo as T* —> 
whichever path is chosen in approaching T* = 0. Even the latter is an extension to the quasi-homogeneous case of 
the standard case [1,67]. See also [68] for further conditions. 



APPENDIX C: GIBBS-DUHEM EQUATION REVISITED 

Let us consider a quasi-homogeneous function g(x 1 , . . . , x n ) of degree q ^ and of weights (a\, . . . , a p , 0, . . . , 0) 
where 1 < p < n is an integer and the last n — p weights are all zero. It satisfies 

Dg = qg, (CI) 

where 



If one defines 

' dx 



dx i 

i<p 



S Vi = l,...,n, (C3) 



and 

= 9i dxi = d 9, (C4) 

i 

one finds 

g=^i D uj g . (C5) 

This formal manipulation is useful in order to derive the Gibbs-Duhem equation. In fact, one has 

dg = d (i i D LUg) = i (-i D duj g + L D u g ) , (C6) 

where u g = q u> g holds. Moreover, du> g = because dg — ui g by construction. The Gibbs-Duhem equation is 
shown to be equivalent to 

--i D du Q = 0. (C7) 
q 

One can easily see that 

duj g = dg- L A dx 1 (C8) 



and that 

i D duj g = Y dgi) dx 1 - ^ d 9i (*r> dx 1 ) 
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= ^2(D gi ) dx l - dgi {oti x*) 

i i<p 

= ^ (l - a i) 9i dx % - ^ d 9i (<*« x 1 ) + ^ 19idx\ (C9) 

i<p i<p p<i<n 

As a consequence, (C7) becomes 

£ - (*T M d ((^) 1 " 9/ai ft) - E 5, d*' = 0. (CIO) 

i<p ^ p<i<n 

If none independent variable of weight zero appears, then one finds 

J2 - (xY a ' d ([x') 1 ^ ft ) = 0. (Cll) 

Let us consider the inverse problem, where one assigns n would-be intensive variables g t and the quasi-homogeneous 
1-form of degree ui g = gi dx l such that Ljj uj g = q u> g (one has to require that deg(gi) — q — ai, where the weights 
of the independent variables x % are the same as above). Then one can also define a function g = ioUg/q- From 
dg = (—id duj g )/q + uj g follows that w g is closed if in dio g = 0, i.e., if the Gibbs-Duhem equation is satisfied. 

We limit ourselves to point out that more general Gibbs-Duhem equations are obtained in the generic symmetric 
case [66]. The generalization of (83) is 

(i x SQrev) d\og ( i ) -i x d SQ rev + (L x - 1) SQ rev = 0. (C12) 



APPENDIX D: QUASI-HOMOGENEITY (SCALING) CANNOT BE GENERALIZED 

We show that, if /(x 1 , . . . , x n ) is a C 1 function on a open connected set which satisfies the following identity: 

f( 9l (X) x\..., 9n (X) x n ) = g(X) f(x\ ...,x n ) (Dl) 

and the functions g(X),gi(\) for i = 1, ...,n are positive definite and invertible, then /(x 1 , . . . , x n ) is necessarily 
quasi-homogeneous. The above transformation is meant to be obtained by a "generalized similarity transformation" 
which carries x 1 , . . . , x n into gi(X) x 1 , . . . , g n (X) x n . Let us consider u = g(X) and the inverse A = <7 _1 (w). Define 

hi{u)= gi {g-\u)) Vi = l,...,n. (D2) 

It is also useful to define 

u = e s . (D3) 

Then we obtain 

/(/M(e s ) x\ . . . , h n {e s ) x n ) = e s f{x\ x n ). (D4) 

Note that, for s = 0, one finds 

/(Ml) x\ . . . , Ml) x n ) = f{x\ x n ). (D5) 

This a priori does not imply that hi(l) = 1 Vz = 1, . . . ,n, but the one-parameter s € IR+ "generalized similarity 
transformation" 7^ : (x 1 , . • • , x") — > (hi(e s ) x 1 , . . . , h n (e s ) x") is consistently defined if 7^, =0 is the identity, i.e. if 
/i,(l) = l Vi = l,...,n. 
We introduce the auxiliary function 

F(s; x 1 , . . . , x") = e- s /(^i(e s ) x 1 , . . . , ^„(e s ) x"). (D6) 

It is such that d s F — 0, i.e. 
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F = J2 e s K(e°) x* (g) (h^) x 1 ,..., K(e s ) x n ), (D7) 
where h'^u) = dhi/du. For s = one obtains 

/(x\...,x")=]T tiJX)* (gW,...,x"). 



(D8) 



Thus, /(x 1 , . . . , x 11 ) is a quasi-homogeneous function of degree one and weights (h[(l), . . . , h' n (l)) = (a\, . . . , a n ). As 
a consequence, one has 

/ (e Ql s x\ ...,e » *x") = e s f(x\ x n ), (D9) 

i.e. 

h,(e s ) = e a > s V* = l,...,n. (D10) 
By resorting the original dependence on A, one finds 

f((g(X)rx\ . . . , ( 5 (A)) a "x") = , 9 (A) f(x\ . . .,x n ), (Dll) 

i.e. 

ffi(A) = ( 5 (A)r Vi = l,...,n. (D12) 

If some <7j = 1 above, then one simply finds that the variable x 1 is quasi-homogeneous of degree zero, i.e., it has weight 
zero. If instead g(X) = 1, then / is quasi-homogeneous of degree zero. The proof is straightforward. Let Ao be such 
that (?j(Ao) = 1 Vi = 1, . . . , n. We have that 

K=0 = f>KA)x* (§P)(gi(X)x\...,g n (X)x n ). (D13) 

By setting A = A , one finds 

E ^(Ao)^ (^)(x\...,x"), (D14) 
i.e., / is quasi-homogeneous of degree zero and weights (g[(Xo), ■ ■ ■ ,g' n (^o))- 

APPENDIX E: SCALING AND ASYMPTOTICS 

Let us assume that f(x) is a function of n variables collectively indicated with x. Let the domain of / be invariant 
under quasi- homogeneous transformations. Let {T\} be a one-parameter quasi- homogeneous transformation. If we 
require that there exist a positive continuous function p(X) and a continuous function g(x) ^ such that g is an 
asymptotic of / in the following sense: 

x^L^x) f{Txx) = a{xl (E1) 

where the limit is assumed to exist. Then: 

a) the function p(X) is a regularly varying function, i.e. it satisfies for all a > 

^lW sC{a) ' (E2) 

where C(a)C(b) = C(a b), i.e. C(a) = a 1 for some 7 <G IR; 7 is also defined the order of the regularly varying function 
P(A); 

b) the function satisfies the "homogeneity" relation 
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g{% x) = u? g(x), (E3) 

i.e. g(x) is quasi-homogeneous of degree 7. 
proof: for each real s > it holds 

A^^** (E4) 



and, because of T\ s = T\ T s , 



As a consequence, 



xToo pj\j /(Ta (Ts X)) = 9{Ts (E5) 



lim J£L JK>*\ = -M- = nm . (E6) 

A^oo p(A s) /(T A (T s x)) g{% x) x^oo p(X s) 



This is possible only if the following conditions 

A^oo p(A) 



fen (E7) 



and 

ff (T s x) = C{s) g(x) (E8) 
are both satisfied. Moreover, it is easy to see that it holds 

C{st) = C{s)C{t), (E9) 
which implies that there exists a real number 7 such that C(s) = s 1 . Thus, one obtains 

g{T s x) = s-> g(x), (E10) 

which concludes the proof. □ 
Cf. also Ref. [55]. 
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